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Abstract 

, We present two loop corrections to photon self energy at finite temperature in real 

(N . time formalism. An expression for renormalized coupling constant has been derived, 

for the first time, in a form that is relevant for all temperature ranges of interest in 
QED, specifically for temperatures T ~ m, where m is electron mass. The latter 
range of temperature is of specific interest from the point of view of the early universe 
' and stellar systems. Such calculations have acquired more significance recently due to 

the possibility of producing e^e~ plasmas in laboratory. We use the calculations to 
determine the dynamically generated mass of photon, Debye screening length, plasma 
frequency up to order as well as the electromagnetic properties of a medium at 
high temperature. 
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> ; 1 Introduction 



Finite temperature field theory (FTFT) wa.s formulated [1-4] a few decades 
. back, for various applications in condensed matter physics, astrophysics and 

cosmology. Research work during the last couple of decades has led to an in- 
creased development and application of methods of using the statistical distribu- 
tion functions in non-relativistic situations [5] as well as relativistic regimes [6]. 
FTFT is relevant to address issues in phase transitions and particle-antiparticle 
symmetry breaking after the Big Bang [7,8] and to trace the thermal history 
of early universe. The temperatures equivalent to electron mass are 
' particularly significant from the point of view of primordial nucleosynthesis and 

QED type corrections are highly important for calculating physically measure- 
able values of QED parameters and their effect on abundance of elements in the 
early universe. 

Finite temperature effects are now extensively used not only in particle 
physics but also to study relativistic plasmas. One of the possible extensions of 
this work is to include the density effect at high temperatures. These extremely 
involved calculations will provide a way to understand super dense electrody- 
namics' systems created in particle colliders and super dense hot stars. Such 
effects are imperative in dense mult i- particle systems in astrophysics and cosmol- 
ogy. These calculations also help to move ahead in understanding phase transi- 
tion in quantum electrodynamics of highly dense systems. More recently, possi- 
bilities to create ultra-relativistic electron positron plasmas with high-intensity 
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lasers has generated renewed interest in finite temperature QED to study pro- 
cesses of pair-production [9-12] in QED plasmas. Lasers pulses hitting a thin 
gold foil can heat up the electrons in the foil up to several MeV's leading to 
pair creation [13-19]. The latest laser facilities plan to attain colossal intensi- 
ties of lO^^W/cm^ which mark a substantial improvement over the ones highest 
recorded (lO^^W/cm^), at Hercules facility in Michigan [20]. 

Quantum Field Theory (QFT) assumes particles to be analogous to excita- 
tions of a harmonically oscillating field which permeates space-time. Modern 
interpretation of vacuum describes it as an absence of particles, but not devoid 
of energy and fields. Vacuum can be treated as a hypothetical bath of virtual 
hot particles which can virtually mediate interactions between real particles for 
un-measureable short intervals of time. Virtual electron-positron pairs couple 
to photons through loops and give photon a temperature dependent finite size, 
charge and mass. Particles propagating in vacuum are assumed to be those 
with interactions switched off. The exact state of all the background particles is 
unknown when these particles propagate through a medium since they contin- 
ually fluctuate between different configurations influencing particle dispersion. 
Thus properties of the system become different from a situation in which all 
particles are assumed to be completely independent of each other, behaving as 
freely propagating bare particles. 

Finite temperature background is particularly interesting as it offers more 
realistic calculation and a possibility of measuring effects of the polarized vac- 
uum by employing thermodynamic quantities. When dealing with sufficiently 
hot environments, the particles are considered to propagate in background heat 
bath at energies around thresholds for particle antiparticle pair production so 
that statistical effects due to temperature need to be appropriately taken into 
account. These effects arise due to continuous exchanges of mediating particles 
during physical interactions that take place in a heat bath containing hot par- 
ticles and antiparticles. Thermal background effects are incorporated through 
radiative corrections and net statistical effects enter the theory through fermion 
and boson distributions. There is multitude of ways to describe physics at finite 
temperature. In real-time formalism, manifest covariance is restored through 
m'* = (1,0,0,0), the four component velocity of background heat bath. Finite 
temperature calculations also provide a guideline to estimate density correc- 
tions at higher order loops through chemical potential effects from background 
plasma. 

Quantum electrodynamics (QED) is useful to study the conditions for exis- 
tence of relativistic electron-positron plasmas. Self energies of particles in high 
temperature medium relevant in QED, acquire temperature corrections through 
virtual exchanges with real particles, in accordance with perturbative nature of 
the theory. Mass less photons gain an effective mass due to medium effects while 
propagating through a system comprising of a cloud of electrons and positrons 
at finite temperature, and vice versa. This gives rise to different terms in pertur- 
bation of states describing the system. Many of these terms can be reproduced 
by modifying particle propagators and from the poles of the propagators one 
obtains modified dispersion relation. Radiatively generated mass shift acts as 
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a kincmatical cut-off in production rate of particles in the heat bath. Gauge 
bosons acquire a dynamically generated mass, at single and higher loop levels, 
due to plasma screening effect. One loop corrections have been studied in detail 
in real time formulation [21-24]. Medium properties such as electric permittiv- 
ity and magnetic susceptibility or permeability also get modified by the thermal 
background effects. 

The problem of rcnormalization in finite temperature theories is similar to 
that at zero temperature. The temperature acts as a regularization parameter 
for ultraviolet divergences. Infrared divergences introduced in finite tempera- 
ture framework are also appropriately removable in particle decay processes, for 
example, via bremsstrahlung emission and absorption effects [22,23]. Renormal- 
ization of QED in this scheme has been already checked in detail at one loop 
level, using real time formulation, for all possible ranges of temperatures and 
chemical potentials [21-29]. The electron self-energy was calculated in detail, for 
all temperature ranges of relevance in QED with renormalized mass and wave 
function determined, at two-loop level [34-37] in this scheme of calculations. 
Method for re-summation over hard thermal loops (HTL) is commonly used to 
determine the higher loop corrections at finite temperature [30-32]. As far as 
rcnormalization is concerned, as discussed in Ref. [33], HTL corrections are not 
necessary. We aim here to calculate photon self energy in a framework valid for 
all temperatures relevant in QED. 

Wc derive an expression for renormalized QED coupling constant in a back- 
ground where electrons and positrons acquire temperature contribution from 
the heat bath along with photons. These calculations have been done in a man- 
ner that the limit of temperature T ^ m, around the threshold for creation of 
e^e~ pairs is inclusive. We now present calculations for photon self energy in a 
general temperature framework, for the first time, for rcnormalization of charge 
in this general framework oiT ^ m so that the calculations are valid for all the 
temperature ranges relevant in QED while the expressions in low temperature 
limit, T « m , presented in Ref. [38], are reproducible from these results. 

Section 2 is based on calculations of second order in a self energy of photons 
that contribute in finite temperature background. Charge renormalization con- 
stant is evaluated in section 3. Expressions for electromagnetic properties in a 
thermal medium are obtained in section 4. Finally section 5 contains discussion 
of the results. 

2 Two Loop Vacuum Polarization 

In real time formalism, temperature dependent (hot) and temperature indepen- 
dent (cold) terms are simply combined, at one loop level, because the propaga- 
tors include them separately as additive terms. At higher-loop level, the loop 
integrals involve a combination of temperature independent (cold) and temper- 
ature dependent (hot) terms which appear due to the overlapping propagator 
terms. Therefore, the calculations are much more involved. 
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Vacuum polarization tensor of photon at the two-loop level contributes to 
the second order hot corrections to cliarge renormalization constant of QED 
at finite temperature. This contribution basically comes from self mass in Fig. 
1(a) and vertex type electron loop corrections inside the vacuum polarization 
tensor in Fig. 1(b). The expression for two loop self energy of photon from self 
energy type correction in Fig. 1(a) is 

(1) 

with k + l = q. On substitution of the particle propagators at finite temperature, 

this gives 

n^,(p,r) = n«,(p,r = o) 



(27r)4y (27r)4 ^"''2ElP[{j> - kf - m?][q^ - w?] 
^ 5'{ko-Ek) ^ 6{ko-Ek) , 6{q^ - m^jn^ (9) 



2 (ko+Ek)' Z2[(p_fc)2 _^2j[;t2 _^2]2 

5[{q - k)^}nB{q ~ k) 



[fc^ — m2]2[(p — k)^ — mPjlq-^ — w?] 

5{{j) - fc)^ - vn?}nF{p - k) 5{{p - fc)^ - vn?}nF{p - k) 

^{q- fc)2[fc2 - m2]2[g2 _ ^2] + _ ^)2[fc2 _ ^2]2[^2 _ ^2] 

2j 5{q^ - m^}nF{q)5{{p - kf - m^jupip - k) 
x6{q - k) nsiq - k)[ + -7j—r^rv]^rE^ 



2 (fco + Ek) ' 2El 

.(5{(p - k)^ - vn?}nF{p - k)5{q^ - m^}nF{q) 

^ (<Z-fc)2[fc2-m2] 

S{{q - fc)^}nB((g - k))S{ip - kf - m^jn^ (p - k) 
[g2 — to2] 



5{{q - kf}nB{q - k)6{q' - m^}nF(g) 
^ (fc2-TO2)[(p-fc)2-m2] J-^J' 



where 
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A/-;, = 8[{3m^ -e -2k.q){ki,{p-k), + K{p-k)^-g^,k.{p-k)} 
+ - m^) {{q - k)^{j) - k), + {q- k),{p - k)^ 
-9t^u{q - k).{p - k)} + 2g^^m'^{m'^ - k.q + k% (3) 

Vertex type correction to two loop self energy of photon in Fig. 1(b) can be 
written as 



(4) 



xf/ I 

[{k -P?- m2] [{q -pY- \q^ - m?] 

5{{q-kf]nB{q-k) _ 5{k'^ - m^)nF{k) 

[fc2 - m2] {q - kf ^ 

1 

[q — kY[{k — pY — rm?] [q^ — m?] 

5{q^ -nY)nF{q) S{{q - pf - nYjupiq - p) . 
[fc2-m2] ^ [g2-rn2] ^ 

S{{k — p)2 - m'^jupik - p) 



+47r2{ 



2r5{{q - kf]nB{q - k)5{q^ - nY}nF{q) 



[q — kY[{q — pY — m?] [k^ — m?] [{q^ — m? 



[{k — pY — 



+ 



5(fc2 - rrY)nF{k) ^ 5{{q - pY — m'^juFiq — p)-. 

[{q — pY — m-2] (/c2 _ jyj2^ 

S{{q - kYjriBiq - k)6{k'^ - rrY)nF{k) 



[q^ — w?] 

r5{{k — pY — m^}nF{k — p) 5{{q— pY — m?}nF{q — p) -, 

[{q — pY ^ [{k — pY — 1^^] [<Y " 

(5(fc2 — m'^)nF{k)5{q^ — m'^}nF{q) 

(5{(g - pY — m^jupiq — p)) S{{k — pY — rrp-^npik, - p) 

\(k — pY — rn"^] [{q — pY ^ '^^l 

5{{q — pY — m^juFiq — p)S{{k — pY — m^juFik — p) 

S{k^ - /)r')//j. (A:) S{q- - ;;/-};) j.-(g) 



where 



^lu = -^[q-{k-p){k^{q-p)„-{q-p)f,K}-q.(q-p){k^{k-p)^ 
+ {k - p)f_ck„} + iq- p)-{k - p){k^q^ + qf_,k^} 

+k.q{{k - p)fj,{q - p)^ + {q - p)^ {k - p)^} - k.{k - p){qn{q - p)v 
+{q - P)iJ-Qi^} + k.{q- p){qy,{k - p)u - {k - p)^,qu} 
+gnt^{q-{k - p)k.{q -p)- k.q{k - p).{q - p) + k.{k - p)q.{q - p)} 
-m^{2fc^(2g - p)„} - [k^{k - p)„ + (fc - + g^uk.{p - k)] 

-[qM - p)v + (9 - p)^Lqu - 9,ii,q-{q - p)\ 

+2{k-p)^,{2q-p),} + m^g^,], (6) 

We find that most of the terms vanish on integrating over hot momenta 
before the cold ones. Therefore, we get respective nonzero contributions to 
these diagrams, only from the terms: 



8{q-kfnB{q-k) 



+47r2<5{g2 _ m''}nF{q)5{{p - kf - m''}nF{p - k)}. (7) 

K^{P, T ^ 0) = STT^e^ / / < <5(9 - kf^B{q - k) 

S{k'^ - m2)nF(fc)(5{(g -p)2 - w2}nF(g-p) 

[(A:-p)2-m2][g2_TO2] ^' ^ ' 

Due to manifest covariance in the theory, physically measurable couplings 
can be evaluated through contraction of vacuum polarization tensor H^^ with 
the metric in Minkowski space g^^" and the bath velocities u^u'^. Following Ref. 
[21], one can write 

n^,(p) = Ul{p)Pi,u + nT(p)Q^., 

where ^ 
and 

nT(p) = -^n4p) + l5''^n^.. 

The coefficients P^i, and Q^,^ are given by 
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and 



P |P| 



If a preferred order of integration (integrating over temperature dependent 
variables before temperature independent ones) is not followed, the loop calcu- 
lations become very cumbersome and the results become messed up. In Eqs. 
(7) and (8), integrating over d'^q (the temperature dependent momentum), us- 
ing Feynman parameterization and then integrating over d'^k, after a somewhat 
lengthy calculation, one gets 



1/1 p2 \ °° (_T\n+r„-nl3E 

.^.-n- = 2^lTHl (i; + 1 - 1^) + 4 \l» _ ^, 

^ (n + 5j m 1 1 

n,r,s — 1 

-r^(4 + ^ In i±^) - 4( In \^flA}]}], (9) 
(r + s) V 1—v v' 1—v 



3 \v 

+48 ^ {(-l)"+'-+i^^[: 



T 

TO 



_ (r — n) (r — n) 

n,r,s — 1 

^ ^ {n + s)\r + s) 
-^{ Ei[-TO/3(r + 5)]}]}], (10) 

with u=Mandi = i- 7- In(^). 

We need to remove the temperature enhanced ultraviolet divergences. The 
counter terms from Fig. (Ic) provide the contribution necessary to cancel the 
divergences appearing in Fig. (la) and Fig. (lb). Adding the results from all 
QED graphs in Figs, (la) - (Ic), it has been explicitly checked that all the hot 
and cold divergences cancel at the two loop level. Therefore, finite terms from 
Eqs. (9) and (10) are used to obtain longitudinal and transverse components of 
the two loop vacuum polarization tensor. Including the one loop corrections in 
Ref. [23] also, these components can be now written as 
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16a, ^ \f/-i ^ 1 + V ..ma(mj3) c{mj3) . 
— (l-^)Ul- 2^in^)( - ^) 

4 ^ 72 
u2^3V 27712^/^ (n + s)(n- 

n,r,s — 

+ V (-l)«+'-e-«^^{-^[-^(4+-lni±^) 
-2m^{ (r + s)/3Ei[-m/3(r + s)]}]}], 



and 



iit{p,t) = -i[ni(p,r)-.9^"^n^,(p,T)] 

8a., 1 1 1 + 11 , ,TOa(TO/3) c(mB). 

= v[^^ + (^-^^^"T^>(^ ^) 



_ (n + s) (r + s) 



--— — In 2m2 - + 6 

[r + sjw 1 — f m 

-(r + s)/3 Ei{-TO/3(r + s)}]}]. 



where 



a(m/3) = ln(l + 6-™^), 

OO 

6(m/3) = ^(-irEi(-nm/3), 



Ia is infrared divergence at finite temperature and Ei is error integral given by 



Ei(-a;) = - 




These expressions are in a form such that all the relevant limits of temperature 
including T ~ m are valid. 



3 Charge Renormalization 

Virtual photons are emitted and reabsorbed by electrons and tend to smear out 
charge distribution and correct the electron mass. The virtual fermion loops 
get polarized in electromagnetic fields, which act to damp the field. Renormal- 
ized electron mass and renormalized wavefunction have been already calculated 
in detail [36]. Vacuum polarization in a medium gives modification to elec- 
tric charge and QED coupling constant. The electric charge couples with the 
medium through vacuum polarization get screened, and picks up thermal correc- 
tions accordingly. We derive an expression for the renormalized charge following 
Ref. [21] and using Eq. (11). The electron charge renormalization constant up 
to the order can therefore be expressed as: 



„ ^ a 8a ,ma(mB) c(md) b(ml3) , , 1 9, 

^3 = 1-^ + ^ - ^^5^^ + + 

67Te 7rm^ p ^ 4 3 



V (-l)^+'-e-^^^-^{- 



w? 6 (n + s) (r + s) 

-12m^[ (r + s)/3Ei{-m/3(r + s)}]}]. (17) 

Using Z3, a corresponding value of QED coupling constant is obtained to be: 



an = a(r = 0){l + ^[— ^ _^ + + _^ )] 



w? 6 {n + s) (r + s) 



g-m/9(s-|-r) 



-12m^[ — (r + s)^Ei{-m^(r + s)}]}]}. (18) 

It can be seen from the above equations that electron charge and hence the 
QED coupling constant gets smaller with an increased order of loops which 
clearly assures the renormalization of electron charge. This modification in 
the coupling constant can be used to determine changes in electromagnetic 
properties of medium. 
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4 Electromagnetic Properties of the Medium 



It is a fact that the electromagnetic properties of media depend on their statis- 
tical properties. Higher order calculation of thermal corrections to self energy of 
photons is much more complicated. However, it is interesting to determine the 
corresponding changes in electromagnetic properties of hot media which arise 
due to radiative emission and absorption of hot particles in the heat bath. It 
can be checked using Ref. [22] that electric permittivity of a medium at the two 
loop level modifies to 



s{p,T) = l-|p|^Hi(p,r) 

IQam? 1 1 + v ma{m(}) c{m(3) 

^^3^ 2m^' {n + s){r-n) 

' ^(n + s)^(r + s) V V \-v 

^-ml3(s+r) 

-2m^{ (r + s)^Ei[-m/3(r + s)}]}], (19) 

whereas magnetic permeability becomes 



8a ,1,^0. 109^;^- 129, ,,,, 



V 



' (n + s)^(s-r) 



[r + s) \ V I — V J \ v 

„-mfi{s+r) 

x{ {r + s)/3Ei[-m/3{r + s)}]}]. (20) 

m 

Equations (19) and (20) show that the radiative emissions and absorptions of 
real hot photons from the heat bath up to order lead to deviations from unity 
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for the values of dielectric constant and magnetic permeability. There arc two 
possible limits for Hl and Ht [22] as p'^ — > and their physical interpretations 
depend on the order of taking the limits po = and |p| — > 0. In Eqs. (11) 
and (12), if we take po = |p| in the rest frame of the heat bath and then limit 
|p| — > is taken, we get 



kI - lim nL(o,|p|,r) 

p| — >o 



16a ma(m/3) c{m(3) m^b{m/3) , 
~^ ^ 2 ^ 



{n + s){r - n) 

s — J- 

4^2 



-m^{n-\-r) 



(n + 5)(r + s) 

n,r,s — ^ 

g-m/3(s+r) 

-2m^{ — (r + s)PEi[-mP{r + s)]}]}, (21) 



4 ~ lim Ul{0,\p\,T) 

IpI— s-O 

-m%{mp) - a-'i V (-l)«+'-_L^{4T^ 

-K ^ ' ^ ^ (n + s) s-r 

n,r,s — 1 

24r / T \ (,-inf){s+r) 

+ 6 mH lOm^ 



(r + s) \ r — s J m 

-{r + s)pEi{-mP{r + s)}]} + ^], (22) 

On the other hand if we set Po = with |p| — > then this result in: 

wl - lim ni(|p|,|p|,T)=0, (23) 
IpI — >o 



^\ - lim nT(|p|,|p|,r) 

Pl 5-0 



8a ,ma(m/3) c{mj3) m?h{ml3) 



-m/3(s-r) / T \ „-m0{s+r) 

+24T— — m + - 12m^' 



{s — r) \ r — s J m 
-{r + s)l3Ei{-mp{r + s)}]}]. (24) 
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The vacuum polarization tensor in — )■ limit provides the dynamically 
generated mass of photon. 

5 Results and Discussion 

Thermal corrections at the two loop level are not as simple as the one loop. 
Temperature dependence to vacuum polarization tensor at the two loop level 
enter as an overlap of photon and fermions hot loop with momenta due to effec- 
tively mutual interactions while propagating through the medium. An unusual 
behavior of hot integrals appears due to the overlap of hot and cold terms. The 
order of integration between the cold and hot loop not only affect the length 
of calculations but also the right order of integration is required to show the 
renormalizability of the theory. With the integration over temperature depen- 
dent variables before temperature independent ones, the loop calculations be- 
come lesser cumbersome and also the cancellation of singularities can be clearly 
seen. In case of wavefunction, the terms with the reverse order of integrations 
have to be included to fully establish the renormalization. The standard reg- 
ularization techniques of vacuum theory like dimensional regularization would 
only be valid in a covariant framework of a Lorentz invariant system. Once the 
hot momenta are integrated out, usual field theoretical techniqiies of Feynman 
parameterization and dimensional regularization can be applied to remove sin- 
gularities arising from cold loops. The results reduce to finite values after order 
by order cancellation of singularities on addition of counter terms of the same 
order. 

At the higher loop level, temperature corrections in hot medium imply con- 
vergence of the perturbativc series. The renormalizcd self energies are used to 
obtain an estimation of dielectric constant and magnetic permeability of a hot 
medium. Both of these quantities deviate from unity even for T<<m whereas 
they do not deviate from the vacuum value with the first order loop corrections 
[23]. This difference of loop correction can be interpreted as at the one loop 
level self-interaction of the photon loops does not contribute whereas at the two 
loop level the hot photon loop develops self-interaction due to the dynamically 
generated mass. This type of effect has been earlier observed for self-mass of 
gluon even at the one- loop level [39] due to self-coupling of gluons. With an 
increase in temperature, photon mass sufficiently affects the behavior of this 
coupling, especially when T > m. This is due to the self mass of photon at 
finite temperature which may lead to nonzero self interaction. 

Further we obtain respective propagation vectors and frequencies when we 
take the limiting values for longitudinal and transverse component of vacuum 
polarization tensor in Eqs. (7) and (8). In particular eqn. (22) nj^^ gives Debye 
screening length of electric force in QED plasma whereas the transverse compo- 
nent of vacuum polarization tensor in Eq. (12) corresponds to the dynamically 
generated mass of photon. QED effective coupling is determined in Eq. (18). 

The framework of calculations adopted here has a form such that the temper- 
atures T ~ m (by taking to/3 ~ 1) are inclusive while the limits of temperature 
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classified as, high temperature {T > m, with e '^^ ignorable as compared to 
^) and low temperature T « m (with fermions contribution negligible) can 
be retrieved from the results obtained in previous sections. 

It is worth-mentioning that all hot corrections give a similar dominant 
dependence, as in case of low temperature [38]. The analysis here is useful 
to evaluate decay rates and scattering cross sections of particles and can be 
generalized for nonzero density that can l)c; used as an effective tool to study 
quark gluon plasma, cores of neutron stars, supernovae, etc. The calculations 
presented here can be extended to the study of electroweak decay rates [40-44] 
at higher loop level. In particiilar such processes will be interesting to study, 
for a better understanding about the medium effects in physical interactions 
associated with black holes and primordial nucleosynthesis. Two loop calcula- 
tions with imaginary time formulation were done; ciarlier for hard lepton pair 
(dilepton) production in the context of quark gluon plasma [45-50]. 
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Figure caption 

Fig. 1. Vacuum polarization diagrams at the second order in a. 
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